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Abstract 

The  air  flow  through  a vane-cascade  swirl  generator  is  examined  both  experimentally  and 
numerically  to  characterize  the  inlet  combustion  air  flow  entering  a reference  spray 
combustion  facility  at  NIST.  A three-dimensional  model  is  used  to  simulate  the 
aerodynamics  in  the  12-vane  cascade  swirl  generator  that  imparts  the  desired  degree  of 
angular  momentum  to  the  air  in  the  annulus  leading  into  the  reactor.  A numerical 
simulation  using  the  Renormalization  Group  method  (RNG)  k-e  turbulence  model  results 
in  a velocity  profile  consistent  with  experimental  measurements,  and  correctly  predicts  a 
recirculation  zone  that  is  experimentally  observed  at  the  exit  of  the  annular  passage.  The 
standard  k-e  turbulence  model  does  not  compare  as  well  with  the  experimental  data  and 
fails  to  predict  the  recirculation  zone  at  the  exit.  This  work  is  part  of  a larger  project  at 
NIST  in  which  benchmark  data  is  collected  for  input  and  validation  of  multiphase 
combustion  models,  and  the  results  presented  provide  a well-characterized  inlet  condition 
for  the  spray  combustion  reactor*.  The  good  agreement  between  the  experimental  data 
and  the  simulation  with  the  RNG  k-e  turbulence  model  provides  further  validation  for  this 
model  in  confined,  annular  flows. 


Keywords : Swirl,  Turbulence,  Fluid  Mechanics,  Numerical  Analysis,  Simulation,  CFD, 
Pitot  Probe 


1.  Introduction 

The  design  and  optimization  of  multiphase  thermal  oxidation  systems  in  the  power 
generation,  waste  incineration,  and  chemical  process  industries  are  relying  increasingly 
on  computational  models  and  simulations  to  provide  relevant  process  information  in  a 
cost-effective  manner.  System  performance  is  dependent  on  the  liquid  atomization, 
chemistry,  aerodynamic  design,  and  the  degree  of  liquid/air  mixing  within  the  reactor.  In 
general,  there  is  a need  for  experimental  data  detailing  the  characteristics  of  the  droplet 
field  and  flame  structure,  and  an  understanding  of  their  interrelationship  with  the  system 
operating  conditions,  for  the  development  and  validation  of  advanced  computational 
models.  This  paper  is  part  of  a larger  effort  at  the  National  Institute  of  Standards  and 
Technology  (NIST)  to  provide  benchmark  data  that  can  be  used  for  input  and  validation 
of  multiphase  combustion  models  and  submodels. 


+ Electronic  files  of  the  data  presented  in  this  report  are  available.  Contact  Cary  Presser  at 
cpresser@nist.  gov. 
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The  experimental  facility  consists  of  an  enclosed  spray  combustion  reactor  and  is 
shown  schematically  in  Fig.  1.  The  combustion  air  passes  through  a 12-vane  swirl 
cascade,  shown  in  Fig.  1 A,  that  imparts  the  angular  momentum  necessary  to  stabilize  the 
flame,  and  flows  around  the  nozzle  before  entering  the  reactor.  The  vane  angle  and  air 
flow  rate  are  50°  and  0.0158  m3  s'1 *,  respectively,  and  correspond  to  the  baseline  case  of 
the  NIST  database  (Widmann  et  al.,  1999).  For  the  reported  air  flow  rate,  the  Reynolds 
number,  Re  = (Do  - D,)u:  /v,  in  the  annular  region  of  the  swirl  generator  is  equal  to  1.0  x 
104.  Note  that  the  radial  locations  of  the  inner  and  outer  walls  of  the  annular  passage 
downstream  of  the  swirl  vanes  are  shown  in  the  figure  (Rt  = 0.0175  m and  Ro  = 0.0508  m, 
respectively). 


EXHAUST 


Fig.  1 Schematic  of  (A)  the  12-vane  cascade  swirl  generator  and  (B)  the  spray 
combustion  reactor.  Dimensions  are  given  in  millimeters. 

In  this  paper,  the  airflow  through  the  vane  cascade  swirl  generator  shown  in  Fig.  1 A 
was  simulated  using  FLUENT+  computational  fluid  dynamics  software  (FLUENT  Inc., 
1998),  and  results  obtained  using  two  turbulence  models  are  presented  and  compared  to 
experimental  data.  The  two  turbulence  models,  the  standard  k-e  model  (Launder  and 
Spalding,  1972)  and  the  Renormalized  Group  theory  (RNG)  k-e  model  (Yakhot  and 

1 Certain  commercial  equipment,  materials,  or  software  are  identified  in  this  publication  to  specify 

adequately  the  experimental  procedure.  Such  identification  does  not  imply  recommendation  or 

endorsement  by  the  National  Institute  of  Standards  and  Technology,  nor  does  it  imply  that  the  materials  or 

equipment  are  necessarily  the  best  available  for  this  purpose. 
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Orszag,  1986),  are  two  equation  models  in  which  scalar  quantities  (the  turbulent  kinetic 
energy,  k,  and  its  dissipation  rate,  £)  are  used  to  describe  the  production,  diffusion,  and 
dissipation  of  turbulence. 

The  standard  k-e  model  is  a semi-empirical  turbulence  model  based  upon  an 
isotropic  eddy-viscosity  hypothesis.  It  is  widely  used  in  industrial  flow  and  heat  transfer 
simulations  due  to  its  robustness,  economy,  and  reasonable  accuracy  (Shyy  et  al.,  1997). 
The  RNG  k-e  model  also  belongs  to  the  k-e  family  of  turbulence  models;  however,  unlike 
the  standard  k-e  model,  the  RNG  k-e  model  was  derived  using  a statistical  technique 
called  renormalization  group  methods.  The  model  equations  are  similar  to  the  standard  k- 
e model,  but  the  statistical  derivation  results  in  different  values  for  the  various  constants 
in  the  equations. 

Smith  and  Reynolds  (1992),  however,  reported  problems  with  the  specific  values  of 
the  constants  in  the  RNG  k-e  model.  In  response,  Yakhot  and  coworkers  reformulated 
the  earlier  derivation  of  the  differential  equation  describing  the  transport  of  £ (Yakhot  et 
al.,  1992;  Yakhot  and  Smith,  1992).  With  this  change,  the  RNG  k-e  turbulence  model 
has  shown  improvement  over  the  standard  k-e  model  when  applied  to  many  industrial 
flows  (e.g.,  Papageorgakis  and  Assanis,  1999;  Yin  et  al.,  1996;  Lien  and  Leschziner, 
1994;  Yakhot  et  al.,  1992).  Particularly  noteworthy  is  the  rate-of-strain  term  in  the 
transport  equation  for  £ that  has  been  reported  to  result  in  improved  predictions  of  flow 
fields  with  high  strain  rates.  In  particular,  flows  in  curved  geometries,  stagnation  flows, 
separated  flows,  and  swirling  flows  are  situations  in  which  the  RNG  k-e  model  has  been 
reported  to  be  more  accurate  than  the  standard  k-e  turbulence  model.  Benim  (1990) 
compared  the  performance  of  these  two  turbulence  models  in  a swirling  combustor  and 
found  that  the  RNG  k-e  model  resulted  in  predictions  consistent  with  experiment,  while 
the  standard  k-e  model  compared  poorly  with  the  experimental  data.  In  this  study,  we 
apply  both  of  these  models  to  simulations  of  flow  through  a vane  cascade  swirl  generator, 
including  the  swirling  flow  through  the  annulus  following  the  vanes.  Jaw  and  Chen 
(1998)  recently  reviewed  second-order  closure  turbulence  models,  and  the  reader  is 
referred  there  for  a more  detailed  discussion  of  these  and  other  turbulence  models. 

2.  Numerical  Methodology 

The  numerical  formulation  is  for  isothermal,  turbulent  airflow  through  the  12-vane 
cascade  swirl  generator  shown  in  Fig.  1A.  The  relevant  conservation  equations  include 
continuity,  the  Reynold's  averaged  Navier-Stokes  equations,  and  an  appropriate 
turbulence  model.  The  Reynolds  averaged  Navier-Stokes  equations  are  generated  from 
the  instantaneous  Navier-Stokes  equations  using  the  following  transformations  (Bird  et 
al.,  1960): 

u.,  = u,  + u\  , p = p + p',  and  p = p + p'.  (1) 

Here,  the  overbar  and  prime  indicate  a time-averaged  quantity  and  an  instantaneous 
fluctuation,  respectively.  Dropping  the  overbar  for  convenience,  the  resulting  expression 
for  the  momentum  equation  in  cartesian  coordinates  is 
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The  last  term  in  Eq.  (2)  is  the  derivative  of  the  Reynolds  stresses,  - pu'u'j , and  represents 
the  effect  of  turbulence  on  the  momentum  balance.  The  Reynolds  stresses  represent 
additional  unknowns,  and  a set  of  constitutive  equations  is  required  to  close  these 
equations.  The  turbulence  models  provide  the  additional  equations  necessary  to  close  the 
transport  equations. 

The  Reynolds  stresses  in  Eq.  (2)  are  computed  using  the  Boussinesq  hypothesis 
(Hinze,  1975), 
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For  the  standard  k-e  turbulence  model,  the  scalar  quantities  k and  £ are  computed  from 
the  following  transport  equations: 
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The  generation  of  k due  to  turbulent  stresses,  Gk,  and  the  generation  of  k due  to  buoyancy, 
Gb,  are  given  by 
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respectively.  Here,  Pr,  is  the  turbulent  Prandtl  number  (Kreith,  1958),  and  /3  is  the 
coefficient  of  thermal  expansion.  For  these  simulations,  the  flow  is  assumed  to  be 
isothermal  and  Eq.  (8)  is  neglected.  The  values  of  the  constants  in  Eqs.  (4)  - (6)  have 
been  determined  experimentally  to  be  C/£  = 1.44,  C2c  = 1.92,  Cp  = 0.09,  Gk  = 1.0,  oe  = 1.3, 
and  Pr,  = 0.85  (Launder  and  Spalding,  1972). 

The  numerical  simulations  were  generated  using  a segregated,  implicit  solver. 
Integrating  the  transient  problem  to  steady  state  was  found  to  be  computationally  less 
expensive  than  solving  the  time-independent  transport  equations,  and  this  method  was 
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used  for  the  results  presented  here.  The  coupled  equations  were  solved  with  first  order 
accuracy  in  time  and  second  order  accuracy  in  momentum,  continuity,  and  turbulence 
parameters.  The  pressure  and  velocity  were  coupled  using  the  PISO  algorithm  (Issa, 
1986)  with  neighbor  and  skewness  correction,  and  standard  wall  functions  (Launder  and 
Spalding,  1 974)  were  used  for  the  near- wall  treatment. 

A uniform  velocity  profile  with  a magnitude  of  2.2  m s'1  was  used  for  the  inlet 
condition,  and  an  inlet  turbulence  intensity  of  10  % was  assumed  with  a characteristic 
length  of  Ro  - R;  = 0.0333  m.  The  predictions  at  the  outlet  of  the  domain  were  not  found 
to  be  sensitive  to  the  inlet  turbulence  intensity.  For  the  conditions  considered,  the 
Reynolds  number  based  upon  the  mean  axial  velocity  is  1.0  x 104  in  the  annular  regions 
of  the  domain.  Sheen  et  al.  (1997)  reported  the  transition  from  laminar  to  turbulent  flow 
to  occur  at  Re  = 1 600;  therefore,  we  expect  fully  developed  turbulent  flow  here.  At  the 
outlet,  the  radial  velocity  was  assumed  to  be  negligible  and  the  radial  equilibrium 
pressure  distribution  was  calculated  by 


dp  _ Pue 
dr  r 

Also  at  the  outlet,  the  turbulence  intensity  and  characteristic  length  used  for  the  inlet 
condition  were  assumed  in  the  event  of  backflow  into  the  domain,  such  as  in  a 
recirculation  zone. 

The  airflow  through  the  12-vane  cascade  swirl  generator  was  simulated  using 
FLUENT  computational  fluid  dynamics  software  (FLUENT  Inc.,  1998).  A three- 
dimensional  model  is  required  for  this  geometry;  however,  due  to  symmetry  it  is  only 
necessary  to  simulate  a 30°  portion  of  the  vane  cascade.  An  unstructured  grid  was  used 
for  the  simulations,  and  the  surface  mesh  is  shown  in  Fig.  2A.  Note  that  the  rotationally 
periodic  symmetry  planes  at  0 = 0°  and  9 = 30°  are  not  shown  in  the  figure.  The  mesh 
was  constructed  so  that  the  grid  resolution  gradually  increased  from  the  inlet  to  the  vanes, 
and  then  remained  high  throughout  the  remainder  of  the  domain.  The  number  of  cells  in 
the  mesh  was  systematically  increased  until  the  solution  was  determined  to  be  grid- 
independent,  and  this  strategy  of  gradually  increasing  the  resolution  from  the  inlet  plane 
to  the  vanes  was  used  for  all  of  the  grids.  The  predictions  presented  here  correspond  to 
results  obtained  from  a grid  with  approximately  277,000  cells.  A top  view  of  the  vanes  is 
shown  in  Fig.  2B,  and  the  30°  section  that  was  modeled  is  depicted.  Note  that  the  vane 
angle,  a,  is  50°  and  shown  in  the  figure. 

3.  Experimental 

The  experimental  facility  consists  of  the  enclosed  spray  combustion  reactor  shown  in 
Fig.  1.  In  this  facility,  the  combustion  air  passes  through  a 12-vane  swirl  cascade,  shown 
in  Fig.  1A,  that  imparts  the  angular  momentum  necessary  to  stabilize  the  flame,  and 
flows  around  the  nozzle  before  entering  the  reactor.  The  reactor  consists  of  a stainless 
steel  enclosure  approximately  1.2  m in  height  and  0.8  m in  diameter,  and  is  shown  in  Fig. 
IB.  The  vane  angle  and  air  flow  rate  are  50°  and  0.0158  m3  s’1,  respectively,  and 
correspond  to  the  baseline  case  of  the  NIST  database  (Widmann  et  al.,  1999).  For  these 
experiments,  the  uncertainties  in  these  measurements  are  ± 1°  and  ± 0.0005  m3  s_1  for  the 
vane  angle  and  air  flow  rate,  respectively.  The  air  enters  the  reactor  at  ambient 
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temperature  and  pressure.  Additional  details  of  the  facility  and  the  benchmark  database 
program  at  NIST  are  available  elsewhere  (Widmann  et  al.,  1999). 


Fig.  2 Schematic  of  (A)  the  surface  mesh  used  to  generate  the  unstructured  grid  for  the 
numerical  simulation,  and  (B)  a top  view  of  the  1 2-vane  cascade  and  the  modeled 
30°  section.  Note  that  the  periodic  boundaries  in  (A)  have  been  omitted  for 
clarity. 


The  magnitude  of  the  air  velocity  was  measured  at  the  outlet  of  the  annulus 
following  the  vanes  as  a function  of  radial  position  to  validate  the  simulations.  A five- 
hole  pitot  probe  with  a hemispherical  head  was  used  for  the  velocity  measurements.  The 
probe  has  a cylindrical  tip  with  a length  of  26  mm  and  a diameter  of  1.7  mm.  Following 
the  26  mm  tip,  the  diameter  of  the  probe  increases  smoothly  from  1.7  mm  to  3.2  mm, 
which  is  the  diameter  of  the  main  body  of  the  probe.  The  probe  has  one  hole  on  the  head 
of  the  probe  tip  for  total  pressure  measurements,  and  four  holes  located  1 3 mm  from  the 
head  and  90°  apart  on  the  side  of  the  probe  for  static  pressure  measurements.  The  four 
velocity  measurements  corresponding  to  the  four  different  static  pressure  holes  were  used 
to  align  the  probe  in  the  direction  of  the  flow,  and  the  average  of  the  four  measurements 
was  reported  as  the  mean  velocity.  The  differential  pressures  were  measured  using 
Autotran  pressure  transducers  (Model  750D-031)  with  measurement  ranges  of  0 - 24.9  Pa 
and  ± 12.5  Pa. 
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To  calibrate  the  pitot  probe,  velocity  measurements  were  obtained  in  a low  velocity 
wind  tunnel  at  NIST.  Data  are  presented  in  Fig.  3 for  four  values  of  the  probe  angle, 
defined  as  the  angle  between  the  direction  of  the  mean  velocity  vector  and  the  probe,  and 
the  error  bars  correspond  to  the  standard  deviation  of  the  four  velocity  measurements.  A 
probe  angle  of  0°  indicates  that  the  probe  was  aligned  with  the  flow.  The  ordinate 
corresponds  to  the  velocity  magnitude  calculated  using  Bernoulli's  equation  (Bird  et  al., 
1960), 


“» =^lPo-P»L  (10) 

where  po  is  the  measured  total  pressure  and  pn  is  the  measured  static  pressure  at  the  nlh 
hole.  The  measurements  obtained  using  Eq.  (10)  are  in  very  good  agreement  with  the 
known  velocity  in  the  wind  tunnel  for  probe  angles  less  than  10°,  and  are  in  fair 
agreement  for  larger  probe  angles.  The  velocity  measurements  reported  in  this  paper 
were  not  computed  using  Eq.  (10),  but  were  determined  from  the  calibration  in  Fig.  3. 


NIST  WIND  TUNNEL  VELOCITY,  m s'1 
Fig.  3 Pitot  probe  calibration  at  four  different  probe  angles. 

To  obtain  accurate  results,  the  pitot  probe  must  be  reasonably  aligned  with  the  mean 
direction  of  the  flow.  However,  Fig.  3 reveals  that  some  deviation  in  the  alignment  is 
acceptable  and  will  not  impact  the  measurements.  For  example,  probe  angles  of  0°  and  8° 
result  in  velocity  measurements  that  agree  within  2 % for  u > 3 m s’1  in  the  calibration 
data.  To  facilitate  aligning  the  pitot  probe  in  the  flow  and  quantifying  the  uncertainties  in 
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the  measurements,  a non-dimensional  alignment  parameter,  'P,  was  used  where  *F  is 
defined  as 


Here,  u is  the  average  of  the  four  velocity  measurements.  The  differences  {uj  - wj)  and 
(w?  - uf)  in  Eq.  (11)  correspond  to  the  velocities  computed  from  the  static  pressure  holes 
that  are  180°  apart  on  the  side  of  the  pitot  probe.  The  alignment  parameter  facilitates 
aligning  the  probe  in  the  flow  by  utilizing  the  symmetry  of  the  probe,  and  should  not  be 
confused  with  a "fractional  error."  However,  'P  was  used  to  estimate  the  uncertainties  in 
the  measurements  by  comparing  data  obtained  in  the  NIST  wind  tunnel  with  data 
obtained  in  the  swirling  flow.  Figure  4 presents  the  calculated  alignment  parameter  from 
the  calibration  data.  Note  that  although  'P  increases  with  increasing  probe  angle  as 
anticipated,  it  is  below  0.07  for  probe  angles  as  large  as  15°.  For  probe  angles  of  0°  and 
8°,  *P  is  below  0.015  for  u > 1.5  m s'1. 


NIST  WIND  TUNNEL  VELOCITY,  m S’1 

Fig.  4 Alignment  parameter,  XP,  calculated  from  the  pitot  probe  calibration  in  Fig.  3. 


4.  Results  and  Discussion 

RNG  k-e  vs.  standard  k-e  turbulence  model 

The  combustion  air  enters  at  the  bottom  of  the  swirl  generator  with  a uniform 
velocity  profile  and  flows  through  an  annular  section  approximately  0.178  m long  with 
the  same  radial  dimensions  as  the  exit  annulus  {Rt  = 0.0175  m and  Ro  = 0.0508  m).  The 
flow  is  then  directed  radially  outward  as  shown  in  Fig.  IB.  The  flow  bends  upward  and 
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then  returns  toward  the  center  of  the  swirl  generator  as  it  passes  through  the  vanes.  The 
swirling  flow  is  directed  upward,  passes  through  an  exit  annulus  approximately  0.165  m 
long  and  enters  the  reactor.  The  fuel  nozzle  is  introduced  into  the  reactor  through  the 
inner  pipe  of  the  annulus;  therefore,  the  swirling  combustion  air  flows  around  the  nozzle 
as  it  enters  the  reactor. 

The  predictions  of  the  magnitude  of  the  air  velocity  through  the  vane-cascade  swirl 
generator  using  the  standard  k-s  and  the  RNG  k-s  turbulence  models  are  presented  in  Fig. 
5.  The  standard  k-s  model  predicts  a maximum  velocity  in  the  annulus  of  ~ 4.8  m s'1, 
while  the  RNG  k-s  model  predicts  a maximum  velocity  of  ~ 5.2  m s'1.  Both  models 
predict  that  the  air  rapidly  accelerates  as  it  enters  the  annulus  following  the  swirl  vanes; 
however,  the  RNG  k-s  model  results  in  larger  velocity  gradients,  with  the  high  velocity 
regions  of  the  flow  being  more  localized.  This  behavior  is  also  evident  in  the  lower 
portion  of  the  domain  beneath  the  swirl  vanes.  The  stagnation  zone  in  the  comer 
corresponding  to  the  first  bend  as  the  flow  leaves  the  lower  annulus  is  larger  for  the  RNG 
k-s  model,  and  all  of  the  recirculation  zones  in  the  domain  have  higher  total  velocities.  In 
comparison,  the  standard  k-s  model  results  in  flow  field  predictions  in  which  the  velocity 
field  is  smoother,  resulting  in  lower  peak  velocities  and  lower  velocity  gradients. 


s 0.0 


STANDARD  k-s  MODEL 


RNG  k-s  MODEL 


Fig.  5 Contour  plots  of  the  predicted  magnitude  of  the  velocity  vector  for  the  standard  k- 
s and  RNG  k- s turbulence  models.  The  contour  plots  correspond  to  the  periodic 
boundaries  at  0~  0°  and  0=  30°,  and  the  scale  ranges  from  0 m s'1  to  5.2  m s'1  for 
both  plots. 
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Contour  plots  showing  the  predicted  turbulence  intensity  in  the  annulus  following 
the  swirl  vanes  for  both  models  are  shown  in  Fig.  6.  Note  that  the  two  plots  have 
different  scales,  with  turbulence  intensities  in  the  range  0.2  - 0.8  and  0 - 0.6  for  the 
standard  k-s  model  and  the  RNG  k-s  model,  respectively.  The  two  turbulence  models 
result  in  significantly  different  predictions  for  the  turbulence  intensity  in  this  swirling 
flow  field.  Both  models  predict  the  generation  of  turbulence  at  the  comer  where  the 
annulus  intersects  with  the  horizontal  section  of  the  swirl  generator  (lower  right  comer  of 
the  contour  plots  in  Fig.  6).  However,  vortex  shedding  is  noticeable  in  the  prediction  of 
the  RNG  k-s  model,  while  the  standard  k-s  model  results  in  significantly  more  dispersion 
of  turbulence.  The  standard  k- s model  predicts  a central  core  of  turbulence  in  the  annulus 
that  decays  with  axial  position.  In  contrast,  the  RNG  k-s  model  predicts  significantly  less 
turbulence  in  the  central  region  of  the  annulus,  with  regions  of  high  turbulence  being 
primarily  in  the  near- wall  region  of  the  inner  wall  and  the  vortex  shedding  near  the  comer 
of  the  outer  wall.  Further  downstream  within  the  annular  passage,  the  flow  separates 
from  the  inner  wall  due  to  the  adverse  pressure  gradient  that  results  from  the  swirling 
flow,  and  which  in  turn  produces  increased  turbulence. 
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STANDARD  k-s  MODEL  RNG  k-s  MODEL 

Fig.  6 Contour  plots  of  the  predicted  turbulence  intensity  in  the  annular  passage 
following  the  swirl  vanes  for  both  turbulence  models.  The  contour  plots 
correspond  to  the  periodic  boundaries  at  6 = 0°  and  0 = 30°,  and  the  scale  ranges 
from  0.2  - 0.8  for  the  standard  k-s  model  and  0 - 0.6  for  the  RNG  k-s  model. 
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Comparison  of  predictions  and  experiment 

Figure  7 presents  values  of  the  alignment  parameter  computed  from  experimental 
data  obtained  at  the  outlet  of  the  annular  region  following  the  swirl  vanes.  The  data 
correspond  to  various  values  of  the  probe  angle  ranging  from  45°  to  60°  measured  from 
vertical;  however,  the  values  of  ¥ at  a particular  radial  coordinate  are  fairly  consistent. 
Near  the  inner  wall,  the  alignment  parameter  becomes  very  large.  It  was  not  possible  to 
align  the  probe  in  the  direction  of  the  flow  in  this  region  due  to  the  presence  of  a 
recirculation  zone;  therefore,  velocity  measurements  were  not  made  in  this  region.  The 
combustion  air  flows  mainly  near  the  outer  wall,  however,  and  the  inability  to  measure 
the  velocity  in  the  recirculation  zone  is  not  considered  significant . 


RADIUS,  m 

Fig.  7 Variation  of  the  alignment  parameter,  'F,  with  radial  position  computed  from  the 
experimental  data  collected  at  the  outlet  of  the  annulus. 

The  predicted  radial  profiles  of  the  velocity  magnitude  at  the  outlet  of  the  annulus 
are  presented  in  Fig.  8.  Note  that  the  location  of  the  inner  and  outer  walls  are  r = 0.0175 
m and  r = 0.0508  m,  respectively.  Experimental  data  obtained  with  the  pitot  probe  are 
also  shown  for  comparison.  The  RNG  k-e  model  compares  well  with  the  data,  while  the 
standard  k-e  model  fails  to  predict  the  recirculation  zone  observed  near  the  inner  wall  of 
the  annulus.  The  recirculation  zone  is  more  apparent  in  Fig.  9 where  the  axial  and 
tangential  components  of  velocity  computed  with  the  RNG  k-e  turbulence  model  are 
presented.  The  radial  velocity  is  negligible  and  is  not  shown.  Near  the  inner  wall  the 
tangential  velocity  is  essentially  zero,  while  the  axial  velocity  is  negative.  It  is  this 


11 


recirculation  zone  that  stabilizes  the  flame  by  transporting  fuel  vapors  from  the 
evaporating  spray  toward  the  reaction  zone.  Thermal  energy  is  also  transported  by  the 
recirculation  zone  in  the  form  of  combustion  products,  thus  contributing  to  flame 
stability.  It  is  therefore  imperative  that  predictive  techniques  capture  this  behavior  if  they 
are  to  be  used  successfully  for  the  design  and  optimization  of  these  processes. 


Fig.  8 A comparison  of  the  total  velocity  profile  at  the  annulus  outlet  predicted  from  the 
numerical  simulation  using  the  standard  and  RNG  k-E  turbulence  models  and  the 
experimental  data. 

The  error  bars  in  Fig.  8 correspond  to  expanded  standard  uncertainties  with  a 
coverage  factor  of  2 (~  95  % confidence  interval  if  the  measurements  are  assumed  to  be 
normally  distributed  about  the  true  mean).  The  procedure  of  Taylor  and  Kuyatt  (1994) 
was  used  in  computing  the  uncertainty  in  the  experimental  measurement.  Uncertainties 
are  classified  as  Type  A or  Type  B,  and  it  is  assumed  that  the  combined  variance  is  the 
sum  of  the  component  variances.  The  Type  A variance  is  a straightforward  calculation 
based  upon  standard  statistical  practices  applied  to  repeated  measurements.  The  Type  B 
uncertainties,  which  can  not  be  quantified  from  replicated  measurements,  were  estimated 
using  the  alignment  parameters  calculated  for  the  data  obtained  in  the  NIST  wind  tunnel 
and  at  the  annulus  outlet.  Variations  in  T*  with  position  and  velocity  were  considered, 
and  data  with  unacceptably  high  values  of  were  discarded.  For  this  reason,  data  close 
to  the  recirculation  zone  have  not  been  reported. 

The  radial  profiles  of  the  turbulent  kinetic  energy,  k,  and  its  dissipation  rate,  £,  at  the 
outlet  of  the  annulus  are  presented  in  Fig.  10.  These  results  were  obtained  with  the  RNG 
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k-e  turbulence  model.  The  turbulent  kinetic  energy  is  relatively  low  in  the  recirculation 
zone,  but  increases  abruptly  in  the  region  of  high  velocity  gradients  at  r = 0.023  m.  Note 
that  the  dissipation  rate,  £,  is  given  on  a log  scale  and  varies  by  six  orders  of  magnitude. 
The  turbulence  intensity  varies  from  ~ 35  % to  ~ 45  % over  the  majority  of  the  annulus 
outlet,  with  lower  values  in  the  recirculation  zone.  Recall  that  the  outlet  of  the  simulation 
domain  corresponds  to  the  combustion  air  inlet  for  the  spray  combustion  reactor  (see  Fig. 
1);  therefore,  the  results  presented  in  Figs.  9 and  10  provide  detailed  reactor  inlet 
conditions  for  the  velocity  components  and  turbulence  parameters,  k and  £,  when 
simulating  the  spray  combustion  process  within  the  reactor. 


RADIUS,  m 

Fig.  9 The  axial  and  tangential  velocity  profiles  at  the  annulus  outlet  predicted  from  the 
numerical  simulation  using  the  RNG  k-£  turbulence  model. 


Swirl  number,  S 

The  degree  of  swirl  present  in  the  combustion  air  entering  a burner  or  furnace  has  a 
strong  effect  on  the  structure  and  stability  of  the  flame,  and  it  is  therefore  an  important 
parameter  in  the  design  and  optimization  of  such  systems.  Despite  the  importance  of 
characterizing  the  swirling  flow,  several  obstacles  prevent  the  reliable  prediction  of 
highly  swirling  flow  fields.  As  discussed  above,  one  difficulty  encountered  when 
designing  these  systems  is  the  questionable  accuracy  of  current  turbulence  models  for 
highly  strained  flows.  This  uncertainty  results  in  flow  field  predictions  that  are  suspect 
until  validated  with  experimental  data.  In  addition,  geometrical  and  empirical 
correlations  available  in  the  literature  are  highly  geometry  dependent  and  often  involve 
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simplifying  assumptions  (e.g.,  inviscid  flow)  that  lead  to  large  uncertainties  in  the 
predictions. 


Fig.  10  Turbulent  kinetic  energy,  k,  and  the  turbulence  dissipation  rate,  e,  at  the  annulus 
outlet  predicted  from  the  numerical  simulation  using  the  RNG  k-£  turbulence 
model. 


Beer  and  Chigier  (1974)  present  a correlation  for  predicting  the  swirl  number,  S, 
obtained  using  vane-cascade  swirlers,  such  as  the  one  used  here.  The  swirl  number  is  a 
non-dimensional  parameter  that  characterizes  the  degree  of  swirl  present  in  the  flow.  It  is 
defined  as  (Gupta  et  al.,  1984) 


(12) 


where  L is  a characteristic  length,  typically  chosen  to  be  the  exit  radius  of  the  burner. 
The  terms  Ge  and  Gz  are  the  axial  flux  of  angular  momentum  and  the  axial  flux  of  axial 
momentum,  respectively,  and  are  given  by 


and 
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Here  uz  and  ug  are  the  axial  and  tangential  velocity  components,  respectively.  For  the  12- 
vane  swirl  generator  investigated  in  the  present  study,  the  characteristic  length  used  was 
radial  coordinate  of  the  outer  wall  of  the  annulus  (see  Fig.  1 A),  L = Ro  = 50.8  mm. 

The  correlation  presented  by  Beer  and  Chigier  is  based  upon  geometrical 
considerations  and  involves  a number  of  simplifications.  Assumptions  made  in  the 
derivation  of  the  correlation  include  inviscid  flow,  a uniform  axial  velocity  distribution, 
and  negligible  angular  momentum  losses  through  the  annular  region  following  the  vanes. 
Note  that  the  pressure  term  in  the  denominator  of  Eq.  (12)  is  neglected.  With  these 
constraints,  the  swirl  intensity  can  be  predicted  from 


where 
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and 


tan(a') 

(1  -X)[l  + tan(a')tan(7r/«v)]  ’ 


(16) 


X = ^ . (17) 

2 nRv  cos  (a1 ) 

Here,  B is  the  vane  height,  A is  the  vane  thickness,  nv  is  the  number  of  vanes,  a'  is  the 
vane  angle,  and  Rv  is  the  distance  from  the  burner  axis  to  the  vanes.  A prime  is  used  to 
denote  the  vane  angle  defined  by  Beer  and  Chigier,  cf,  and  to  differentiate  it  from  the 
vane  angle,  a,  defined  in  Fig.  2B.  Figure  11  shows  the  relation  between  a'  and  a.  The 
angle  a is  measured  from  the  center  of  the  vane  as  shown  in  Fig.  2B,  while  ot  is 
measured  from  the  edge  of  the  vane.  In  the  above  equations,  ct  is  the  ratio  of  the  average 
tangential  and  radial  velocity  components  exiting  the  vanes,  and  X is  a blockage  factor 
that  accounts  for  the  thickness  of  the  vanes. 


Fig.  11  Relationship  between  the  vane  angle,  a,  and  the  vane  angle  as  defined  by  Beer 
and  Chigier,  c£. 
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Using  Eqs.  (15)  - (17),  the  calculated  swirl  number  for  the  swirl  generator 
considered  here  (a  = 50°,  a'  = 70.7°)  is  S = 1.0.  The  value  of  the  swirl  number 
determined  from  Eqs.  (12)  - (14)  and  the  results  of  the  numerical  simulation  (RNG  k-e 
turbulence  model)  is  S = 0.49.  Each  of  the  assumptions  presented  above  will  lead  to  a 
geometrical  correlation  that  overpredicts  the  swirl  intensity.  In  addition,  a comparison 
between  experimental  data  and  the  correlation,  Eqs.  (15)  - (17),  shows  that  the  swirl 
intensity  is  overpredicted  by  the  correlation,  and  the  discrepancy  between  the  theory  and 
data  increases  with  increasing  values  of  a (see  Fig.  5.4  in  Beer  and  Chigier,  1974).  The 
experimental  data  presented  do  not  contain  values  of  a larger  than  1.5,  while  the  value  of 
a computed  for  the  swirl  generator  used  in  this  study  is  2.1;  thus,  one  would  expect  the 
actual  value  of  S to  be  considerably  less  than  that  predicted  by  the  correlation,  as  is  the 
case  here. 
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Fig.  12  Variation  of  the  swirl  number  in  the  annulus  with  axial  position.  The  swirl 
number  was  calculated  from  the  simulation  using  the  RNG  k-e  turbulence 
model. 

Equations  (15)  - (17)  predict  the  swirl  intensity  immediately  following  the  vanes, 
and  therefore  neglect  the  decay  of  angular  momentum  that  occurs  in  the  annular  passage 
downstream  of  the  swirl  vanes.  The  effect  of  neglecting  angular  momentum  losses  can 
be  explored  using  the  numerical  simulation  by  evaluating  S at  various  axial  locations  in 
the  annulus.  Figure  12  presents  the  computed  values  of  the  swirl  number  in  the  passage. 
The  swirl  intensity  of  the  air  entering  the  annulus  immediately  downstream  of  the  vanes 
is  approximately  S = 0.54,  and  it  is  reduced  to  S = 0.49  by  the  exit  of  the  annulus.  While 
the  swirl  number  is  reduced  due  to  the  angular  momentum  losses  in  the  annular  passage, 
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it  does  not  appear  to  be  the  dominant  reason  that  the  geometrical  correlation  overpredicts 
the  exiting  swirl  intensity  for  this  configuration. 

The  swirl  number,  S,  has  been  shown  to  be  an  important  similarity  variable  for 
geometrically  similar  swirl  generators  (Beer  and  Chigier,  1974;  Gupta  et  al.,  1984). 
However,  it  should  be  used  with  caution  when  comparing  the  swirl  intensity  produced  by 
different  types  of  swirl  generators.  There  is  some  ambiguity  due  to  the  choice  of  a 
characteristic  length  in  Eq.  (12),  and  S is  often  calculated  neglecting  the  pressure  term  in 
Eq.  (14),  which  can  result  in  much  higher  values  of  S.  It  is  therefore  prudent  to  use  care 
when  designing  or  optimizing  processes  based  upon  swirl  number  correlations. 

5.  Conclusions 

The  air  flow  through  a vane-cascade  swirl  generator  was  investigated  both 
experimentally  and  computationally.  The  experiments  reveal  a recirculation  region  near 
the  inner  wall  at  the  outlet  of  the  annulus,  and  this  region  is  predicted  using  the  RNG  k-e 
turbulence  model.  The  standard  k-e  turbulence  model  is  unable  to  predict  the 
recirculation  region  for  this  flow  field,  and  results  in  overly  dispersive  velocity  and 
turbulence  predictions.  Previous  investigators  have  reported  superior  performance  for  the 
RNG  k-e  model  over  the  standard  k-e  model  for  highly  strained  flows,  and  the  results  of 
this  investigation  provide  further  validation  for  this  turbulence  model  for  confined, 
swirling  flows. 

The  simulation  incorporating  the  RNG  k-e  model  results  in  a swirl  number,  S,  of 
0.49  exiting  the  annular  region  following  the  swirl  vanes.  A geometrical  correlation  from 
the  literature  overpredicts  the  swirl  intensity  by  a factor  of  two,  and  the  likely  causes  have 
been  discussed.  In  addition,  the  inlet  condition  for  the  combustion  air  corresponding  to 
the  baseline  case  of  NIST's  benchmark  database  for  the  validation  of  multiphase 
combustion  models  has  been  characterized.  Readers  interested  in  obtaining  electronic 
files  of  the  data  presented  herein  should  contact  Cary  Presser  (e-mail: 
cpresser@nist.  gov). 
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